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ABSTRACT 

The annihilation spectrum of pairs with one-dimensional thermal distributions in the 
presence of a strong magnetic field is calculated. Numerical computations of the spectrum 
are performed for mildly relativistic temperatures and for different angles of emission with 
respect to field lines. Teragauss magnetic fields are assumed so that conditions are typical 
of gamma-ray burst and pulsar environments. The spectra at each viewing angle reveal 
asymmetric line profiles that are signatures of the magnetic broadening and redshifting 
of the line: these asymmetries are more prominent for small viewing angles. Thermal 
Doppler broadening tends to dominate in the right wing of the line and obscures the 
magnetic broadening more at high temperatures and smaller viewing angles. This angular 
dependence of the line asymmetry may prove a valuable observational diagnostic tool. For 
low temperatures and magnetic field strengths useful analytic expressions are presented 
for the line width, and also for the annihilation spectrum at zero viewing angle. The 
results presented find application in gamma-ray burst and pulsar models, and may prove 
very helpful in deducing field strengths and temperatures of the emission regions of these 
objects from line observations made by Compton/GRO and future missions. 

INTRODUCTION 

There is substantial observational motivation for detailed calculations of the annihi- 
lation spectrum of thermal pairs in the presence of a strong magnetic field. If the features 
observed around 400 keV in gamma-ray burst (GRB) spectra (see Mazets et al . , 1981) are 
interpreted as red-shifted 511 keV lines (for an extensive list of observations see Liang, 
1986), then pair annihilation is taking place near a neutron star surface and will thus be 
influenced by a strong magnetic field. Further, recent observations (see Agrinier et al., 
1990) by the FIGARO II experiment of an emission feature at about 400 keV in the spec- 
trum of the Crab pulsar arc suggestive of a red-shifted annihilation line in this neutron 
star source. A study of the influence that a strong magnetic field has on the annihilation 
spectrum is therefore of significant interest to astrophysicists. 

There are two important effects of the magnetic field on two-photon pair annihilation: 
(i) because of the very high cyclotron emission rates, the vast majority of pairs will annihi- 
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late from the ground Landau state with a one-dimensional (ID) momentum distribution, 
and (ii) since transverse momentum is not conserved in a strong magnetic field, there is an 
intrinsic broadening of the two-photon annihilation line due to an increased availability of 
phase space for this interaction. These effects become significant when the magnetic field 
exceeds about a few percent of the quantum critical field B c = m 2 c 3 / (eh) = 4.413 x 10 13 
Gauss, which corresponds to conditions in neutron star environments. To a first approxima- 
tion, it is reasonable to assume that the annihilating pairs have a ID thermal distribution. 
The effects of the field on the shape of the thermal annihilation spectrum are potentially 
the most useful observational diagnostic tool: the total rate of annihilation is not altered 
much by the presence of the field when B 4 x 10 12 Gauss (Wunner, 1979). 

Due to the intrinsic geometrical asymmetry imposed by the field, both the thermal 
and magnetic contributions to the width of the annihilation line will be very angle depen- 
dent. Intuitively, at viewing angles 0\ ~ 0° to the field, thermal Doppler broadening is 
expected to be a maximum while magnetic broadening is a minimum. For viewing angles 
~ 90° to the field this situation is reversed so that the magnetic broadening should be 
more significant in comparison with the Doppler broadening. These expectations are borne 
out by the computational results that are presented here. The prediction that the width of 
observed lines at some angles could be due primarily to the magnetic field is observationally 
significant. Annihilation spectra produced by pairs with a ID thermal distribution have 
been calculated by Kaminker et al. (1990), revealing significantly reduced Doppler broad- 
ening at 6\ ~ 90° . However, they did not include the effects of magnetic broadening. In an 
earlier paper, Kaminker et al (1987) presented spectra for the annihilation of pairs at rest 
in the presence of a field. Here we present the first calculation of two-photon annihilation 
spectra from a ID thermal pair distribution, including full magnetic broadening effects 
using the QED second-order annihilation cross section of Daugherty and Bussard (1980, 
hereafter DB80). Of particular interest are the angle-dependent asymmetries of the line 
profiles introduced by field broadening, signatures that may be detectable by the BATSE 
instrument aboard the GRO mission and future high-resolution detectors (e.g. TGRS and 
IVAE). 


MAGNETIC TWO-PHOTON ANNIHILATION 

It is convenient to use a simple dimensionless notation for various quantities in rel- 
ativistic magnetized environments, in order to express analytic results more compactly. 
Hereafter, all magnetic fields will be expressed in units of the critical field B c , and the 
parallel pair temperature parameter T will represent the quantity k B T e /(m e c 2 ) , where 
k B is Boltzmann’s constant. All photon and particle energies and momenta will be scaled 
by m e c 2 and m e c , respectively. 

The conservation relations for energy and also for momentum parallel to the magnetic 
field during pair annihilation can be simply written down. Let p_m e c and p+m e c be the 
momenta of the annihilating electron and positron, oppositely directed along the field, and 
£_m e c 2 and e+m e c 2 be their energies. The produced photons have energies uqm e c 2 and 
u> 2 m t c 2 , and propagate at angles 8 X and 8 2 to the field. Using the abbreviated notation 
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Cj = cos Oi , Si = sm 6i , i = 1,2 , the conservation relations are 


c- + 64 - = + ^2 

p_ — p+ = Cc?iCi + C^2 c 2 


( 1 ) 


Alternatively these relations can be expressed in terms of the solutions for the pair 

momenta p± (see DB80). The kinematical requirement that these momenta be real gives 
a restriction on the energy-angle phase space of the emitted photons. 


Ul\ C\ + i^2 c 2 


< 


\j (uq + ^ 2 ) — 4 


(2) 


For the annihilation of pairs with momenta p±m e c along the field, the distribution 
function for the pairs that is needed for the computation of thermal annihilation spectra 
can be obtained from equation (23) of Harding (1986): 


**±(P±) = 2I\i(l/T) exp { *T } ’ 

where T is the dimensionless temperature of the pairs and K x is a modified Bessel func- 
tion. The thermal annihilation spectrum can then be calculated by inserting this into the 
expression in eq. (24) of DBSO, who performed the QED calculation for the annihilation 
cross-section from first principles. Their result is too small by a factor of two (c.f. Wun- 
ner, 1979), and also contains some minor typographical errors. The corrected annihilation 

spectrum taken from DBSO is then 


where 


Q 27 


2\2 

n n _ c Qf A c uj\ 
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uq -f- W2 ^2 
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(4 b) 


Here o f w 1/137 is the fine structure constant, A c = h/(m e c) and n± are t ^e particle 

number densities. The matrix element M can be obtained from the results of DBSO. The 
sum is over the polarization states of the photons. 

In obtaining eq. (4) a change of variables was performed: the energy u, 2 and angle 
e 2 of the unobserved photon were expressed in terms of the momenta p± of the pairs via 
eq (1) This transformation is convenient for both numerical and analytic computations, 
yielding simple integration limits: with it is associated a Jacobian given by u 2 du 2 dc 2 - 
|p_g , -f p + £_ !/(£_£+) dp-dp + . The azimuthal angle cf> 2 of the second photon is the third 

integration variable. In the limit of small 9 X , the dominant contribution to the annihilation 
spectrum occurs when the exponential factor in eq. (4b) is maximized. 


247 


NUMERICAL RESULTS 


The results of numerical computations of the spectrum in eq. (4) are displayed in 
Figures 1-3. The principle magnetic effects are demonstrated in Fig. 1, where the pairs 
are assumed to be at rest, i.e. zero temperature. This corresponds to the case studied 
numerically by DB80 and Kaminker et al. (1987), and indeed the results presented in 
Fig. 1 are in close agreement with these earlier studies. With thermal broadening absent, 
Fig. 1 clearly illustrates the magnetic broadening effect in the left wing of the two-photon 
annihilation line. There are two prominent features to notice. First, the widths of the 
annihilation lines are dependent on the strength of the field and the angle of emission 
relative to field. These magnetic widths can be expressed approximately (for B <C 1 ) as 



si < \/2B , 
si > \/2B , 


( 5 ) 


remembering that si = sinfh . These estimates can be obtained from eq. (4) or from the 
results of Ivammker et al. (1987). Second, observe the appearance of angular-dependent 
kinematic cutoffs, eliminating the right (blue) wing of the line. These occur at 


“ irki ’ ( 6 ) 

as can be determined from eq. (2). Due to the presence of this kinematic cutoff at uq > 1 , 
most of the broadening occurs on the left (red) side of the line at angles other than 90° . 
Note that a summation of the spectra over emission angles removes the evidence of the 
line asymmetry and the kinematic cutoffs, producing a symmetric line profile. 

The inclusion of thermal broadening for the case of finite T is illustrated in Figs. 2- 
3. The basic nature of Doppler broadening is similar to the ID field-free case (Kaminker 
et al . , 1990), and the results show that the thermal contribution to the line broadening is 
greatest for small viewing angles. The first order contribution to the Doppler width is 

(Auq) D = aV2T , (7) 

when T < 1 . The thermal broadening smooths out the line profile, and dominates the 
blue wing of the line, but evidence of the magnetically-broadened left wing remains for 
larger viewing angles. A comparison of eqs. (5) and (7) reveals that magnetic broadening 
will dominate over first-order Doppler broadening when 


/ 0.25 B s\jc\ , si>s/2B, 

\ 0.125 B 2 /c], si<s/2B 


confirming the numerical spectral results that the Doppler mechanism is more effective at 
broadening the line for small viewing angles. 

The calculations performed here are for conditions expected to be typical of gamma- 
ray burst or pulsar environments. It is anticipated that pairs in the ground state are likely 
to be in relatively cool quasi-thermal distributions. Cyclotron resonant scattering can 
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FIGURE 1 - Two-photon annihilation spectra for pairs at rest in the ground state, for 
different magnetic fields B and viewing angles #i from numerical integration of the cross 
section over the energy u >2 and angle 62 of the unobserved photon. The kinematic cutoffs 
occur when | cos 82 \ — 1 (cf. Eq. 6). 



FIGURE 2 - Numerical annihilation spectra at an angle of 0° to the field direction for 
pairs in a ID thermal distribution with parallel temperature T (in units of m e c 2 ) at two 
different field strengths. 
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FIGURE 3 - Numerical annihilation spectra at different viewing angles 0 1 and parallel 
pair temperatures T, showing the effect of combined magnetic and thermal broadening. 

cool pairs in the ground state before they annihilate (Preece and Harding, 1991) to tem- 
peratures T ~ B/A (Lamb et al, 1990). In this case magnetic broadening is dominant for 
tan 2 0i > 1 or at viewing angles 0j > 45° . This angular effect shows significant potential 
as a spectral diagnostic for line observations in astrophysical sources. 
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ANALYTIC SPECTRAL FORMS 

It is useful to develop analytic approximations to the annihilation spectrum. These 
can be simply obtained in the low T , low B limit, and when the viewing angle is small: 
« 0° . The integrations in eq. (4) are then analytically tractable because of the effective 
restriction to the o >2 — cos 82 (or equivalently the p_ — p+ ) phase space. F 01 low T , we set 
p_ = p+ = 0 in M . Further, specializing to Si « 0 simplifies M dramatically (only one 
term of a series remains — see DB80 for details). In this case, A4 — A^ 1 ^ + A^ 2 ^ , with 


A^> = 

W (2) = 


2 

B +u> 1 

2 

B + u> 2 




( 9 ) 


and the polarization vector components e z , e-j- are given in DB80. The int.egiation o\er 
azimuthal angle turns out to be trivial, yielding 


R 27 



2 2 
2B 


£— + £ 4 - — 

T 



(10a) 


yy iioi v 2 2 

a = tb “* E{K (i, r + K ( 2 > r} ■ do*) 

pol 

Since the pair temperature is low, the integrations over p_ and p + can be performed by 
the method of steepest descent. However it must be noted that there are two major regions 
of contribution to these integrals: in the neighbourhood of p_ - -p+ = 0 , corresponding 
to c 2 = —w\C\jbJ 2 and u> 2 = 2 — U\ , and also near c 2 = -1 and u) 2 = l/u>i . The total 
spectrum is, in general, well approximated by the sum of these two contributions. 

The contribution R, from p± « 0 for = 0° occurs only in the left wing of the 
line. Evaluation of eq. (10) is then straightforward, resulting in 


7?i(u)i) ~ 


n+n. 


.ca f 2 A 2 


2a .’ 1 a) 2 — 2aq 4- 2 
~B~ (R + aj ) 2 



2(1 -o>i) 

B 


ui < 1 


( 11 ) 


This contribution is magnetically broadened, and is generally the dominant contribution 
when T & B 2 . It is consistent with the numerical plots, and reproduces the exact cold 
plasma, zero viewing angle approximation derived in eq. (14) of Kaminker et al. (1987). 
The second component B u , for C 2 ~ — 1 , p+ ~ 0 and p_ ~ oq — 1/aq , occurs in both 
wings of the line and is a Doppler broadened contribution. When B ^ \/T , eq. (10) gives 

RuM » n+n-cafXl expj- ^"^ , (12) 


which agrees with the B = 0 , ID thermal result in eq. (30) of Raminker et al. (1990). 
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The results in equations (11) and (12), which can be extended to provide a more 
complicated approximation to the spectrum when T ~ B 2 , indicate that the annihilation 
line is thermally broadened in the right wing ( uq > 1 ), but both thermally and magnet- 
ically broadened in the left wing. This concurs with the numerical results plotted in the 
figures. The widths obtainable from eqs. (11) and (12) agree with those in eqs. (5) and (7), 
and from eq. (8) it can be deduced that when 5 < 1, small temperatures T B 2 /8 are 
necessary for magnetic broadening to be important when 9\ « 0° . Work is in progress to 
extend these analytic results to the regime of larger angles 9\ . 

CONCLUSIONS 

The results presented in this paper show that two-photon annihilation of electrons 
and positrons with ID thermal distributions can be significantly influenced by a strong 
magnetic field. Even for temperatures as high as T = 0.1 , at large viewing angles to the 
field, magnetic broadening dominates for B 10 12 G. More specifically, in gamma-ray 
burst models where the cooling of pairs by cyclotron resonant scattering is important, 
the widths of annihilation features seen at viewing angles 6 > 45° will be determined 
primarily by the magnetic field. At most viewing angles, the magnetic field preferentially 
broadens the red side of the annihilation line. An integration over viewing angles removes 
this asymmetry. Analytic approximations for zero viewing angle are obtained: the spec- 
trum is well- described by a simple sum of magnetically-broadened and Doppler-broadened 
contributions. For low T , the numerical computations are found to agree well with the 
cold annihilation results of Kaminker et al. (1987). The results presented in this paper 
have significant observational implications: the widths and spectral shapes of observed 
annihilation features in gamma-ray burst and pulsar spectra may provide diagnostics of 
both the magnetic field strength and viewing angle in these sources. 
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